Stretched exponential relaxation in a diffusive lattice model. 
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We studied the single dimer dynamics in a lattice diffusive model as a function of particle density 
in the high densification regime. The mean square displacement is found to be subdiffusive both 
in one and two dimensions. The spatial dependence of the self part of the van Hove correlation 
function displays as function of r a single peak and signals a dramatic slow down of the system for 
high density. The self intermediate scattering function is fitted to the Kohlrausch- Williams- Watts 
law. The exponent f3 extracted from the fits is density independent while the relaxation time r 
follows a scaling law with an exponent 2.5. 
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I. INTRODUCTION 

Lattice-gas models of dense fluids have been explored 
Jl|-[| in order to understand the microscopic details of 
their slow dynamics and to compare them with the be- 
havior of other disordered systems like structural glass 
formers, spin glasses and colloids. 

Among the lattice models Random Sequential Adsorp- 
tion (RSA) of particles on a substrate is of particular 
interest. RSA has been invented to mimic many ex- 
perimental situations like chemisorption on crystal sur- 
faces, reactions on polymer chains, adsorption of macro- 
molecules and colloidal particles, monolayer and mul- 
tilayer growth processes (for an an extensive review 
see ||,B)). In particular RSA with added Diffusion 
(RSAD) has been considered a model which captures the 
main features of the compaction dynamics in granular 
matter || . The description of the relaxation processes in 
the RSAD model has been investigated through the be- 
havior of the time density correlation function in a lim- 
ited number of studies and only during the deposition 
process j^]. 

In this paper we explore by computer simulation the 
equilibrium dynamical properties of a lattice gas of 
dimers as a function of density in two dimensions. The 
starting configurations have been produced with a RSAD 
model where by deposition and diffusion of particles we 
reach the desired value for the density. In this way, due 
to the geometrical frustration, as the dimer density in- 
creases trapping configurations that only depend on the 
local environment || are created in the system. The 
dynamical properties are studied by switching the depo- 
sition off. 

In the next section we describe the model and the de- 
tails of our Monte Carlo simulation. In the third section 
we present our results for the single dimer dynamics as 
a function of particle density. In particular we analyze 



the mean square displacement, the Van Hove self corre- 
lation function and its space Fourier transform, the self 
intermediate scattering function. Last section is devoted 
to concluding remarks. 



II. MODEL AND COMPUTATIONAL DETAILS 

In order to investigate the dynamical properties 
of dimers on a square lattice with fixed density we 
have generated the starting configurations through a 
RSA model of dimers with diffusion in two dimensions 
(RSAD2d) §,0 by Monte Carlo (MC) simulations. We 
performed also simulations of the same model in one 
dimension (RSADld) for comparison. We have consid- 
ered a square lattice of N = L x L sites, with periodic 
boundary conditions in both directions. Two kinds of 
dynamical processes can occur in the system, namely 
deposition of a dimer on an empty pair of sites and diffu- 
sion of a deposited dimer on an empty first neighbor site; 
no overlap of dimers is allowed. The lattice is initially 
empty and at each MC step one deposition event (or one 
diffusion event) is selected with probability p (or 1 — p), 
where < p < 1. Dimers can be adsorbed on the lattice 
in two possible orientations (horizontal and vertical) and 
can diffuse along their own axis. In the case of deposi- 
tion, the dimer can occupy a randomly selected pair of 
sites if they are both empty, otherwise the attempt is 
rejected. In the case of diffusion the randomly selected 
horizontal (vertical) particle attempts to move left or 
right (up or down), with equal probability. A diffusion 
event to the right (left) or up (down) is possible only if 
the chosen pair of sites is occupied and the corresponding 
right (left) or up (down) neighbor site is empty, otherwise 
the attempt is rejected. As usual the MC time unit corre- 
sponds to N deposition/diffusion attempts. We checked 
that finite-size effects are negligible for L > 40, so in 
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the following we analyze a square lattice with L = 40. 
We produced 11 starting configurations at densities p = 
0.60, 0.70, 0.75, 0.80, 0.88, 0.94, 0.96, 0.98, 0.99, 0.995, 0.998 
For each configuration production runs at constant den- 
sity have been obtained by switching the deposition 
events off. In this way we are able to follow the dynam- 
ics of the system in an equilibrium situation. In fact, 
the equilibrium distribution for this system, where all al- 
lowed configurations have the same weight, is a constant. 
Therefore the detailed balance condition reduces to two 
requirements: transitions to forbidden states are not al- 
lowed, and transition from an allowed state to another 
has the same probability as the reverse one. Both these 
requirements are satisfied by the adopted algorithm. 

A more complete investigation of the irreversible dy- 
namics of these models and discussion of the role of the 
dimensionality is presented elsewhere || . 



III. SINGLE DIMER DYNAMICS 



In this section we investigate the dynamical behavior 
of this lattice model through the single dimer density 
correlators. We have chosen relatively large values of the 
density [p > 0.60), since we are interested in the diffusive 
dynamics of the system in the high densification regime 
where slow dynamics is more likely to occur. The time 
over which the diffusive motion of particles has been fol- 
lowed reaches in most cases 10 6 MC time units (because 
of the statistical noise in the last decade only t up to 10 5 
is shown in the plots). 



A. Diffusion properties 



We analyze in the following the mean square displace- 
ment (MSD) of the particle, < r 2 (t) >, defined as 
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FIG. 1. Log-log plot of the temporal evolu- 
tion of the MSD for several different densities: 
p = 0.60, 0.70, 0.75, 0.80, 0.88, 0.94, 0.96, 0.98, 0.99, 0.995, 0.998 
from top to bottom. Long dashed lines are the fit to power 
law behaviour in the two time regions from 5 to 120 MC 
timesteps and from 1000 to 60000 MC timesteps. For the 
sake of clarity fits are reported only for even curves. The 
exponents extracted from the fits are reported in table |l| 

After an initial transient, which is more pronounced 
at high densities, the temporal behavior of the MSD can 
be fitted with two power laws. In the first time region 
< r 2 (t) >~ t Sl . In the intermediate time region for 
the lower densities the MSD shows a bent followed by 
a slight increase in the exponent, i.e. < r 2 (t) >~ t 52 
with $2 > Si, where the exponents depend upon density. 
The bent becomes more evident as the density increases. 
Values of <5i and 82 are reported in table 0. 



TABLE I. Values of the exponents extracted from the fit 
to the power law of the MSD plotted in Fig.l. 



Density 


Si 


82 


S2/S1 


0.60 


0.75 


0.80 


1.07 


0.70 


0.72 


0.80 


1.11 


0.75 


0.67 


0.79 


1.18 


0.80 


0.64 


0.79 


1.23 


0.88 


0.60 


0.74 


1.23 


0.94 


0.54 


0.70 


1.30 


0.96 


0.52 


0.67 


1.29 


0.98 


0.51 


0.60 


1.18 


0.99 


0.51 


0.48 


0.94 


0.995 


0.51 


0.38 


0.75 



N part 

< r 2 (t) >= < £ \ ri (t) ri (0)| 2 >, (1) 

part . 1 



where N part = N ■ p/2 is the number density, ri is the 
position of particle i at time t and <> denotes a time 
average. In fig. |l| we show a log-log plot < r 2 (t) > for 
RSAD2d for the different densities investigated. 
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They are less than 1 for all the densities considered: 
in particular 5i and 82 do not exceed respectively 0.75 
and 0.80, which are the values found for the lowest 
density, i.e. p — 0.60. The difference between the 
two exponents is more marked for intermediate densities 
(0.80 < p < 0.96). For the highest density (p = 0.998) 
the MSD shows a plateau for large times. In the time 
region investigated none of the curves shown attains the 
Brownian diffusive regime, 82 = 1- 

These findings can be interpreted in terms of the den- 
sity characterizing the degree of motion of particles like 
the thermal energy in thermal systems. For low den- 
sities single particles can move freely and the diffusive 
dynamics is almost independent of time. Increasing the 
density particles remain trapped in local transient config- 
urations PJul and diffusion slows down. For very high 
density the trapping configurations are frozen and the 
system does not evolve any more on the time scale of 
our simulations. The existence of trapping configura- 
tions and subdiffusive behavior is in general linked to the 
presence of disorder that in the present case can be 
ascribed to the random barriers associated with the local 
environment determined by the geometrical constraints. 
The particle has to overcome these barriers and spend 
a long time before detrapping (see [Q for an analogous 
case in a random walk model). 

These dynamical features are not present in the 
RSADld model, where we have found a MSD with a 
density independent slope for this case. Nevertheless the 
behavior is subdiffusive, with an exponent S = 0.5 con- 
stant over the whole time window investigated. 

B. Density correlator 

For a closer analysis of the motion of the dimer we 
have considered the self part of the van Hove correlation 
function (SVHCF) G s (r,t), defined as 

AW* 

G.(r,t) = -—< <5(r-N(i)-ri(0)|)>, (2) 

ls part i=1 

which gives the probability of finding a particle around a 
distance r at time t, given that the same particle was 
at the origin at time t = 0. The MSD analyzed in 
the previous subsection is the second spatial moment of 
the SVHCF. Nonetheless, in disordered diffusive systems 
that do not have a simple behavior, the SVHCF might 
contain additional information on the microscopic relax- 
ation mechanisms of the system, that could be smeared 
out in the MSD. 

The SVHCF is shown for our system in fig. || for four 
densities at different times. We have calculated the one 
dimensional SVHCF along the horizontal or the vertical 
direction. In fig. || we show the SVHCF averaged over 
the two directions. The self correlation functions shows a 
single peak centered around zero for all the times investi- 
gated. The appearance of more than one peak at a given 



time would have evidenced the existence of more than 
one, well distinct, relaxation mechanism for that time. 




FIG. 2. Spatial dependence of the SVHCF G s (r,t) for dif- 
ferent densities: p = 0.60 (a), p = 0.80 (b), p = 0.88 (c) and 
p = 0.98 (d). The curves in each graph correspond to various 
times: t = 3,15,63,127,511,1023 (in MC time units) from 
top to bottom on the left side of each graph. Distances on 
the horizontal axis are measured in lattice spacing units. 

The single peaked structure confirms the dynamics to 
be mastered only by a single mechanism related to the 
existence of high entropic barriers generated by configu- 
rations in which a particle is locally trapped. It is also 
seen that at low density the SVHCF is localized at short 
distances (typically 1 < r < 2) for short times and tends 
to delocalize for large times, since it is able to explore 
more configurations. At higher densities we observe a 
"clustering" of the various curves for 1 < r < 2.5 (this 
effect is particularly evident in fig. |^(d)): this is the signa- 
ture of the fact that the movement of particles drastically 
slows down at such densities. 



C. Intermediate scattering function 

The space Fourier transform of the SVHCF, i.e. the 
Self Intermediate Scattering Function (SISF) F s (q,t), is 
given by 

AW* 

F s (q,t) = - expItq-MtJ-riCO))], (3) 

where q = (q x ,q y ) is the discrete wave vector (q x , y — 
2nn X: y/L with n x , v an integer between and L — 1). For 
all the densities considered in the simulations F s (q,t) is 
found to decay more slowly the smaller q is. We show in 
fig. H data for q = (2n(L — 1)/L,0) at several different 
densities. This q value is chosen of the order of the in- 
verse of the typical size of the trapping region, which we 
estimate to be of the order of 2 — 3 lattice spacing M. 
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As it can be seen, the entire dynamical behavior at low- 
intermediate densities is observed in the time window of 
our simulations, since the SISF decays to zero for such 
densities. For high values of the density only a part of 
the SISF can be observed, due to the very slow dynamical 
relaxation and the freezing of particles which is also evi- 
dent in the MSD. The short-intermediate time behavior 
of SISF can be fitted with a Kohlrausch- Williams- Watts 
law (KWW) 

F a (q,t) ~ exp[-(tM% Q</3<1. (4) 

The curve resulting from this fitting, also shown in fig. 
reproduces the data obtained from MC simulations quite 
accurately for several decades. 




t 



FIG. 3. Temporal evolution of the SISF F s (q,t) for 
q = (2n(L — 0). The solid curves correspond 

to the result of the simulations for different densities: 
p = 0, 60, 0.70, 0.75, 0.80, 0.88, 0.94, 0.96, 0.98, 0.99, 0.995, 0.998 
The long-dashed curves are fits to the data according to 
eq. (^). In the inset it is shown the inverse relaxation time 
vs. density: the squares are the parameters obtained from 
the fit to eq. (3) for various density, while the solid curve is 
the best fit according to eq. (JE]). 

A deviation of MC results from eq. (Q) is observable 
only in the tails of the correlators. The collapse of cor- 
relators onto a single curve when plotted versus t/r{p) 
demonstrates the existence of a well defined scaling law 
on several decades for this kind of systems. In the range 
of density considered the stretching exponent (3 is essen- 
tially constant, with a value ~ 0.68, the density depen- 
dence being reflected only in the relaxation time r. 

Furthermore we have found that the behavior of r ver- 
sus p follows the scaling law 

T-^A^-pr. (5) 

This fit is shown in the inset of fig. || together with the 
points determined from the parameters of eq. (Q). As it 
can be seen from the figure, eq. (^) applies satisfactorily 
in the entire range of densities considered. Our estimates 
yield A ~ 0.0023, 7 ~ 2.5 and p c ~ 1.0. Given the un- 
certainties in the fit, this value found for p c cannot be 



distinguished from the maximal density reachable in this 
system, which has been found to be s=s 0.998 in numerical 
simulations (see |0,Eol)- This would suggest that there 
is no dynamical transition in this model, i.e. no critical 
density (lower than the maximal density) exists at which 
a structural arrest of the system occurs. It is important 
to stress that the relaxation time obeys the power law 
eq. (^|) , which has the same kind of scaling found in many 
glassy systems and is the signature of a sluggish dynam- 
ical behavior. The existence of a power law behavior for 
t as a function of density and the validity of KWW law 
are novel and interesting points for systems of this kind. 
It would be interesting to investigate in the future the 
t and [3 dependence on the wavevector to check on the 
range of validity of both the KWW law and of the power 
law. 



IV. CONCLUSIONS 



We presented a Monte Carlo study of the single dimer 
density correlators for a randomly diffusive model of 
dimers on a square lattice as function of density. For each 
density the starting configurations have been generated 
through a RSAD2d model. We found that the dynamical 
behavior is related with the system dimensionality and 
that the diffusive dynamics is non Brownian for the time 
regimes investigated. In spite of its simplicity the model 
considered shows a quite complex dynamical behavior, 
and in particular some features resembling that of the 
slow relaxation of structural glass formers are found. 

The most interesting results came from the SISF anal- 
ysis. We have found that this correlator in the q, t space 
has a well defined scaling behavior with density that can 
be described with a high degree of accuracy by the KWW 
function in most of the time region investigated. This 
function is well known to reproduce many features of 
disordered systems, and is connected to fractality in con- 
figurational space Jl3|,|l4| ■ The exponent of the stretched 
exponential function [3 is density independent and the re- 
laxation time follows a power law with an exponent 7 in 
the range of those of structural glass formers . The 
power law predicts a divergence of the relaxation time 
for density around one. From these findings it can be 
consistently derived that for this non thermal system the 
role of thermal energy is played by the density. 

We correspondingly find a subdiffusivc behavior in the 
MSD. The slow down of the system that becomes more 
significant as density is increased can be attributed to 
transient trapping configurations. We point out that 
with respect to lattice-gas models of the kind considered 
in ref. M we did not introduce any explicit constraint to 
induce the trapping of the dimers. In our system in fact 
it naturally arises because of geometrical constraints. 
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